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NOTATION 

Area  of  section 
Coefficient  of  area 
Coefficient  of  I 

o 

Coefficient  of  I 

^  o 

Lift  coefficient 

Thrust  coefficient  for  nonviscous  flow 
Diameter  of  propeller 

Maximum  ordinate  of  mean  line  (camber  line) 

Moment  of  inertia  about  an  axis  parallel  to  the  a>axis  and  passing  through  the  center  of 
gravity 

Moment  of  inertia  about  the  z-axis 

Moment  of  inertia  about  an  axis  parallel  to  the  y-axis  and  passing  through  the  center  of 
gravity 

Moment  of  inertia  about  the  y-axis 
Length  of  blade  section 
Bending  moment  from  torque  per  blade 
Bending  moment  from  thrust  per  blade 

Bending  moment  about  abscissa  with  reference  to  the  center  of  gravity 

Bending  moment  about  ordinate  with  reference  to  the  center  of  gravity 

Pitch  at  any  radius 

Torque  per  blade 

Torque  for  nonviscous  flow 

Radius  of  propeller 

Radius  of  mean  line 

Variable  radius 

Radius  of  any  propeller  blade  section 

Length  of  mean  line 

Maximum  thickness  of  section 

Thrust  per  blade 

Thrust  for  nonviscous  flow 

Speed  of  advance  of  propeller 
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x'  Distance  along  mean  line 

x  Abscissa  of  center  of  gravity  with  reference  to  midpoint  of  section 
ajj  Abscissa  of  nose  with  reference  to  axis  through  center  of  gravity 
*2  Abscissa  of  tail  with  reference  to  axis  through  center  of  gravity 

*3  Abscissa  of  point  of  maximum  thickness  with  reference  to  axis  through  center  of  gravity 

y'/ 2  Ordinate  of  the  section  measured  perpendicular  to  the  mean  line 

y  Ordinate  of  center  of  gravity  with  reference  to  nose»tail  line 

yx  Ordinate  of  nose  with  reference  to  axis  through  the  center  of  gravity 

y2  Ordinate  of  back  with  reference  to  axis  through  the  center  of  gravity 

y3  Ordinate  of  point  of  maximum  thickness  with  reference  to  axis  through  the  center  of 
gravity 

2  Number  of  blades 

Hydrodynamic  pitch  angle 
e  Drag-lift  ratio 
0  Angle  of  arc  of  the  mean  line 
p  Density 

d>  Tan- 1  (— — )  =  Corrected  hydrodynamic  pitch  angle 
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ABSTRACT 


An  approximate  method  of  obtaining  stress  in  a  propeller  blade  from 
simple  beam  theory  is  presented.  The  method  is  designed  to  minimize  the 
work  required  for  calculating  the  geometric  properties  of  the  blade  section. 

An  example  is  given  with  this  method  applied  to  a  modern  thickness  form. 

INTRODUCTION 

The  stress  calculation  usually  applied  to  ship  propellers  was  developed  by  D.W.  Taylor. 
The  basic  assumptions  are  that  the  simple  beam  theory  is  applicable  to  a  propeller  blade  and 
that  thrust  and  torque  depend  linearly  on  the  radius.  Usually  the  maximum  stress  is  determined 
for  only  one  section  of  the  blade  close  to  the  hub.  The  thickness  of  the  sections  is  then 
assumed  to  vary  linearly  with  the  distance  from  the  hub. 

As  to  the  first  assumption,  a  closer  approximation  to  the  stresses  can  be  obtained  if 
the  propeller  blade  is  considered  as  a  cantilever  plate2  rather  than  as  h  simple  beam.  How* 
ever,  the  time  required  for  determining  the  stresses  of  a  cantilever  plate  appears  to  be  ex¬ 
cessive  for  the  present  application.  Therefore,  the  assumption  of  the  simple  beam  theory  is 
retained  in  the  following  consideration.  As  to  the  dependence  of  thrust  and  torque  on  the 
radius,  the  circulation  theory,  developed  since  Taylor’s  work,  makes  possible  the  calculation 
of  thrust  and  torque  distributions  on  a  propeller  blade.  Assuming  these  distributions  to  be 
known,  the  bending  moments  and  stresses  at  any  section  can  be  calculated;  thus  the  minimum 
thickness  of  a  section  at  any  radius  may  be  determined. 

The  problem  then  becomes  one  of  obtaining  the  different  geometric  properties  of  the 
blade  sections.  Muckle3  has  developed  a  simple  method  of  obtaining  the  center  of  gravity 
and  moments  of  inertia  for  ogival  sections.  This  method  enables  the  designer  to  obtain 
minimum  thickness  for  these  sections  wi  thout  numerical  integration  for  the  geometric  proper¬ 
ties  of  the  section.  It  is  important  that  the  blade  thickness  be  kept  at  a  minimum  for  efficien¬ 
cy  and  cavitation  considerations;  therefore,  there  is  a  need  for  a  method  of  obtaining  the 
geometric  properties  of  sections  of  any  shape  quickly  and  accurately.  In  principle,  these 
properties  can  be  derived  by  means  of  numerical  integration;  however,  this  is  time  consuming. 

An  approximate  method  is  developed  in  this  report  for  determining  the  stress  in  a  propel¬ 
ler  blade  and  an  example  is  given  of  its  use.  A  method  for  determining  the  geometric  proper¬ 
ties  of  a  blade  section  (similar  to  that  given  in  Reference  3)  is  given  for  the  TMB  EPH  thick¬ 
ness  form  superimposed  on  several  mean  lines.  By  introducing  correction  factors,  this  has 
been  applied  to  other  thickness  forms  in  use  today,  i.e.,  NACA  16  and  NACA  65A  superim¬ 
posed  on  the  same  mean  lines. 


^References  are  listed  on  page  18. 
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STRESSES  FROM  BEAM  THEORY 


The  stresses  in  a  propeller  blade  are  caused  by  thrust,  torque,  and  centrifugal  force. 
The  most  important  stresses  are  the  bending  stresses  due  to  the  thrust  and  torque.  Bending 
stress  due  to  centrifugal  force  requires  consideration  when  the  blades  are  raked.  The  other 
stresses,  that  is,  shear  caused  by  thrust  and  torque  and  tension  caused  by  centrifugal  force, 
are  usually  negligible. 

For  this  particular  development,  only  the  bending  stresses  caused  by  thrust  and  torque 
will  be  considered  (Figure  1).  For  the  determination  of  stress  caused  by  centrifugal  force 
see  References  1  and  3. 

Assuming  that  the  thrust  and  torque  distribution (d  Tb/dr  and  dQb/dr)  on  the  propeller 
blade  are  known,  the  bending  moment  at  radius  rQ  caused  by  the  thrust  is  given  as: 


dTk 


M 


Tkm  I  (f-'o )-£~df 


and  that  by  the  torque  as: 


where  R  is  the  radius  of  the  propeller  and  r  is  a  variable  radius. 

The  bending  moments  are  resolved  into  two  moments:  M  about  an  axis  parallel  to  the 

* 6 

nose-tail  line  of  the  section  and  the  other  M  perpendicular  to  this  line.  Both  of  these  axes 

xQ  and  yQ  pass  through  the  center  of  gravity  (centroid)  of  the  blade  section  (Figure  2).  When 

the  pitch  angle  <f>  and  the  moments  and 

b 

are  known,  the  bending  moments  about 
the  a?0-  and  yQ- axes  are  given  as  follows: 

cos  <4  +  sin  $ 


M 


v  -*  M  >r  sin  <4  -  Mn  cos  <4 
“o  ii>  Vfc  r 


The  positive  sense  of  each  of  these  moments 

is  indicated  in  Figure  2. 

If  /  and  /  are  the  moments  of  inertia 
o  yo 

about  the  *0-and  yQ-axes,  the  stresses  in  the 
blade  are  given  by  the  following  relationships: 


Figure  1  -  Forces  on  a  Blade  Section 
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Stress  in  leading  edge 


Stress  in  trailing  edge 


y,  M  x  Mv 
*  9  *  7  o 


/. 


I 


?o 


y$  mx 

Stress  on  back  at  point  of  maximum  thickness  ■  - _ 2 

lx 

o 

A  positive  stress  denotes  tension  and  a  negative  stress  denotes  compression. 

As  shown  in  Figure  2,  ®l,  x2>  and  a?3  are  abscissas  of  the  nose,  tail,  and  point  of 
maximum  thickness,  respectively.  Likewise,  y j,  y2,  and  y3  are  ordinates  of  the  nose,  tail, 
and  point  of  maximum  thickness.  It  should  be  noted  that  the  x -  and  y-coordinates  must  be 
used  with  their  proper  signs  in  order  to  obtain  the  correct  signs  for  the  stresses  from  the 
above  formulas. 

• 

QUANTITIES  ARISING  FROM  THE  GEOMETRY  OF  THE  SECTION 

*  Once  the  bending  moments  are  derived,  the  problem  becomes  one  of  obtaining  the  geo¬ 

metric  properties  of  the  blade  sections.  These  properties  are  evaluated  by  solving  the  follow¬ 
ing  integrals  between  their  proper  limits  in  which  the  a-axis  is  the  nose-tail  line  and  the 
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Figure  3  *  Coordinates  of  Blade  Section 

y-axis  is  a  line  perpendicular  to  and  passing  through  the  midpoint  of  the  nose-tail  line  (Fig¬ 
ure  3). 

Area 

A  =  //  dydx 

Horizontal  distance  to  center  of  gravity 

ffx  dydx 
**  A 

Vertical  distance  to  center  of  gravity 

_  ffv  dvdx 

y - J“ 

Moment  of  inertia  about  9-axis 

lxx  -  //y2  dyd* 

Moment  of  inertia  about  y-axis 

lyy  -If*2  dVdx 

These  integrals  are  evalutated  for  several  section  forms. 
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TABLE  1 

Offsets  for  TMB  EPH  Sections 


x7 1 

CM  | 

x'/l 

y  72 
t 

0 

0 

0.500 

0.4946 

0.005 

0.0756 

0.550 

0.4830 

0.010 

0.1064 

0.600 

0.4647 

0.020 

0.1498 

0.650 

0.4399 

0.040 

0.2092 

0.700 

0.4085 

0.070 

0.2717 

0.750 

0.3705 

0.100 

0.3186 

0.800 

0.3260 

0.150 

0.3774 

0.850 

0.2749 

0.200 

0.4204 

0.900 

0.2170 

0.250 

0.4522 

0.930 

0.1778 

0.300 

0.4750 

0.950 

0.1480 

0.350 

0.4902 

0.970 

0.1129 

0.400 

0.4983 

0.990 

0.0642 

0.43613 

0.5000 

1.000 

0 

0.450 

0.4997 

Nose:  0^x7^0.4361302 


x7*  \»  /ym» 
0.4361302/  \  t  / 


=  1 


Body:  0.4361302^x7^0.8722604 


\( 

2  \0.43 61302 


1 


Tail:  0.8722604^x7^1 


1  -  x7 1 
0.3  0  83906 


The  radius  of  the  nose  is  0.5732  t2/7r  the  radius  of  the  tail  0.2027  1 2/J 
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EPH  ■  THICKNESS  FORM,  CIRCULAR  ARC  MEAN  LINE 

The  usual  method  of  approach  is  to  solve  the  preceding  integrals  by  the  use  of 
Simpson’s  Rule.  A  direct  solution  can  be  made  if  the  section  chosen  can  be  expressed 
mathematically.  The  TMB  EPH  thickness  form  (see  Table  1)  is  formed  of  parts  of  an  ellipse, 
a  parabola,  and  a  hyperbola;  therefore,  it  can  be  expressed  mathematically.  Likewise,  a 
circular  arc  mean  line  can  be  represented  by  a  simple  mathematical  expression.  For  this 
reason,  the  foregoing  integrals  wer,e  solved  using  the  EPH  section  in  combination  with  a  cir¬ 
cular  arc  mean  line. 

The  method  of  solution  was  first  to  transform  the  formulas  both  for  the  three  parts  of 
the  EPH  section  and  for  the  geometric  properties  into  polar  coordinates  r  and  9  (Figure  3). 

Thus,  using  polar  coordinates,  the  formulas  are  expressed  as  follows,  where  y'/ 2  is 
the  ordinate  of  the  section  measured  perpendicularly  to  the  mean  line: 

02<*c+y'/2)  0, 

A  f  f  rdrdQ  =  f  y'RedQ  [1] 

(*c-r'/2) 

d2  (Vr'/2)  e2  .  ,2v 

Sin  *  Sin  *4  *  i2  n?r 


/ffc2  sin20  y'dO 

i®.  ,  when  t/l  ^0.21  and  f/l  ^0.05 

•  — a - 


02  (Rg+y'm 

f  f  r[f  cos  9  ~(Re  ~f)]  drdd 

6x(Re-y'm 

V - ~A -  t3] 

A 


l*xm  f  f{Rc+y  /2)r[r  cos  9 /)]*  drdd  [4] 

9X  (Rc  -y'/2) 


o2  ^+y'/2)  e2  . 

Kym5  /  ^sin2  e  drdd  -  /  R3  sin20  y'  (  1  +  i  —  )<J0 

9l(Re-y'/ 2)  9t  '  *Re2' 


2 

■  /  Re3  sin20  y'dd,  when  t/l% 0.21  and  f/\%. 0.05 


[5] 
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ft 


Within  the  range:  0.06387  —  ^0  =  0.5 


ft 


ft 


t 


0.43613 


J/o.18613  +  0. 


ft0  ft202 

12774  — - e- — 

s  a 2 


Within  the  range:  -0.37226  —  ^  0^  0.06387  — 


ft. 


ft 


[6] 


/  ff  0  ft2#2  \ 

y '  =  2.62868  « ^0.37634  +  0.12774  -2 - J 


Within  the  range:  -  0.5  —  <0<  -  0.37226  -5- 

Re~  ~  Rc 


0.61678 


r 


55839  +  1.61678 


ft0  R}92 


a  a-1 


[7] 


[8] 


In  these  equations,  Re  is  the  radius,  f  is  the  maximum  height,  and  a  is  the  length  of 
the  mean  line  and  l  is  the  length  of  the  nose-tail  line. 

The  next  step  involves  substituting  Equations  [6],  [7],  and  [8]  into  Equations  [1],  [2], 
[3],  [4],  and  [5].  It  will  be  noted  that  when  these  substitutions  are  made,  integrands  of  the 
form  sin  n6\Jf(d)  and  cos  "0  \fJT&)  are  obtained.  An  approximate  solution  for  these  integrals 
can  be  obtained  by  expanding  sin  0  and  r~  o  0.  Assuming  sin  0  ■  0  -  03/6  and 
cos  0  =  1  -  02/ 2  +  04/24,  a  very  close  approximation  will  be  obtained  since  the  camber  ratio 
f/l  is  normally  less  than  0.05.  For  this  value  of  f/l ,  sin  0  and  cos  0  are  correct  to  four  places. 

This  integration  is  performed  separately  for  each  of  the  three  parts  of  the  blade  section. 
By  adding  the  results  of  these  separate  integrations,  the  geometric  properties  of  the  whole 
section  are  obtained  as  follows: 

A  =  0.746  at 
r-  0.027  a  +  0.00528 

*c2 

y  =  f  -  0.0256—  +  0.06  ~ 

Rc  Rc 

l„  «  0.04487  a  t3  +  0.08975  iil/  +  0.001626 

XX  Rc  R2 

+  0.746  atf2  -  0.04468  +  0.00188  ?-L 

Rc  R? 
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/  -  0.04469  a3*  -  0.012943 

vy  R? 


For  use  in  the  stress  formula,  the  moments  of  inertia  must  be  taken  about  the  center  of 
gravity  and  xv  *2,  xy  yy  y2,  and  y3  must  be  obtained  from  x  and  f  (Figure  2).  Thus  the 
following  equations  result: 

A  -  0.746  at 

x1  -  0.52  -  0.027*  -  0.00528  a3/R2 


x 


2 


-  I 


»  -0.5 1  +  <0.5 1  +  Re )  sin  MfiMT  * 

R  e 

y.  -  0.0256  ~  ~  0.060  £  -  f 


y2  = 


y3  -  (0.5 1  +  Re)  cos  1  +  f  —  R e  +  y} 

K  C 

L  -  y2  ■=  0.04487  <2  +  0.00391  ^ 

°  L  ftc2 


-  0.00645  1^  +  0.00139  0.00269 


R . 


R- 


a  t 


lyy  -  AS2 


0.04415  -  0.01298  —  -  0.000021 1^- 


R. 


an 


Since  from  geometry 


a 


1 

3 


(4  y/l  +  4  (f/l)2 


small  cambers), 


R  -  fiiMl  i 

c  [  8  f/l 

and  noting  that  t  =  (t/2)  l  and  f  =  ( f/l)l ,  the  equations  can  be  written  in  terms  of  the  section 
length,  thickness  ratio,  and  camber  ratio. 

These  equations  are  cumbersome,  but  their,  use  can  be  simplified  by  limiting  the  thick* 
ness  ratio  to  the  maximum  value  of  0.21  and  the  camber  ratio  to  0.05.  The  equations  are 
then  evaluated  for  a  number  of  different  thickness  and  camber  ratios.  By  plotting  the  resulting 
values,  the  formulas  may  then  be  reduced  to  the  following  forms: 


A  =  (a  t/l  )l2  a  from  Figure  4 


=  (0.47325  -  0.026  f/l)l 

x2  “  *i  *  * 

*3  -  *1  -  0.4361  l  +  0.17  l 
yx  *  -  (0.114  «/i  +  0.79)  (//I)  i 

y2~yl 

y3  «  0.5  t  +  0.984 /+  yl 


QOI 4 
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Figure  5  -  Curve  for  Obtaining  Moment  of  Inertia  (/  )  of  Blade  Section 

x  xr\ 
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Figure  6  -  Curve  for  Obtaining  Moment  of  Inertia  ( ly  )  of  Blade  Section 

The  quantities  obtained  for  the  EPH,  circular  arc  section  have  been  numerically  com¬ 
pared  with  the  respective  quantities  for  the  EPH  thickness  form  combined  with  a  NACA 
a  *=  1.0  and  a  =  0.8  mean  line.  Within  the  limits  of  accuracy  of  integration,  the  different  mean 
lines  had  no  effect  on  the  numerical  results  obtained. 

The  appendix  gives  an  example  of  the  use  of  the  above  equations. 

NACA  16  THICKNESS  FORM,  «  =  0.8  MEAN  LINE 

The  NACA  16  thickness  form4  has  an  elliptical  nose  similar  to  the  EPH  section  but 
with  a  thinner  tail.  Using  this  thickness  form  with  a  NACA  a  =  0.8  mean  line,  the  geometric 
properties  were  obtained  numerically  for  various  thickness  and  camber  ratios.  From  these 
results,  correction  factors  were  determined  for  the  formulas  developed  previously.  Thus,  the 
equations  become: 

A  -  0.886  {a  t/l)  l 2  a  from  Figure  4 
xt  •  (0,4838  -  0.026  f/l)  l 

a?.,  =  -  l 

A  L 

X,  i:t  •  0.8/ 


yx  *-(0.113///  +  0.T82)  {f/l)  l 
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y2-y  i 

y3  =(o.5<  + 

/  =  0.9925  [5  (//l)2  +  0.04487  {t/l  )3  ]  l*  b  from  Figure  5 

o 

l  -  0.946  {c  t/l  )Z4  o  from  Figure  6 
yo 

When  the  quantities  obtained  for  this  section  (NACA  16,  a  «  0.8)  are  compared  numeri¬ 
cally  with  those  obtained  for  NACA  16,  circular  arc  and  a  -  1.0  sections,  it  can  be  shown 
that  the  different  mean  lines  have  no  effect  on  the  numerical  results. 


NACA  65A  THICKNESS  FORM,  a -0.8  (MODIFIED)  MEAN  LINE. 

This  section  is  a  modified  form  of  the  NACA  65  section.  A  mean  line,  NACA  a  »  0.8 
(modified),  was  designed  for  this  form.3  The  NACA  65A  thickness  form  is  desirable  for  ship 
propellers  since  the  nose  of  the  section  is  more  blunt  than  an  ellipse,  and  the  tail  is  thicker 
than  that  of  the  NACA  65  section. 

Using  the  same  method  as  followed  with  the  NACA  16,  a  «  0.8  section,  the  formulas 
for  the  geometric  properties  are: 


A  -  0.903  (a  t/l)  l 2  a  from  Figure  4 


*1  -(0.4467  -0.026  f/l)l 


X 2  “*|  ~l 


*3  “  * 


t  -  0.41  +  0.26  l 
yx  -  -  (0.115  t/l  +  0.798)  {f/l)  l 


V2mV  i 

y3  -  0.5$  +  0.96 /+  yj 

L  -  0.864  [6  {f/l)2  +  0.04487  (t/l)3]l4  b  from  Figure  5 
o 

/  -  0.789  (c  t/l )  lA  c  from  Figure  6 

'  o 


As  with  the  EPH  and  NACA  16  sections,  the  various  mean  lines  have  little  effect  on 
the  numerical  results. 


SUMMARY 

1.  It  has  been  assumed  that  the  simple  beam  theory  can  be  applied  to  a  propeller  blade. 
Certain  approximations  can  also  be  made  to  facilitate  integration  when  determining  the 
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geometric  properties  of  the  section. 

2.  The  formulas  for  the  geometric  properties  of  the  EPH  section  are  obtained  in  simplified 
form  by  limiting  the  thickness  ratio  to  0.21  and  the  camber  ratio  to  0.05.  These  formulas 
apply  to  the  EPH  section  when  used  with  the  circular  arc  mean  line.  These  same  formulas 
may  also  be  used  with  the  mean  lines  NACA  a  =  1.0  and  a  *=  0.8  with  sufficiently  accurate 
results. 

3.  Formulas  for  use  with  the  NACA  16  and  NACA  65A  thickness  forms  are  obtained  by 
applying  correction  factors  to  the  formulas  for  the  EPH  section.  The  resulting  equations  are 
applicable  when  these  sections  are  used  in  conjunction  with  the  circular  arc,  NACA  o  -  1.0, 
NACA  a  «  0.8,  or  NACA  a  -  0.8  (modified)  mean  lines. 

4.  If  desired,  correction  factors  may  be  obtained  for  a  number  of  other  sections  by  a  meth¬ 
od  similar  to  that  used  here. 
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APPENDIX 

SAMPLE  STRESS  CALCULATION 


The  following  calculations  apply  to  a  propeller  blade  with  EPH  sections  and  circular 
arc  mean  lines.  It  will  be  assumed  that  the  distributions  of  the  thrust  coefficient  (dCT)/dx, 
lift  coefficient  cL,  and  hydrodynamic  pitch  angle  p.  are  known  from  propeller  theory.  The 
particulars  for  the  propeller  are  assumed  as  follows: 

Diameter  -  20  ft 
Speed  of  advance  *  50  fps 
Number  of  blades  =  6 
Length  of  section  =  5.0  ft 
Camber  ratio  *  0.013 
Thickness  ratio  «*  0.18 
Pitch  at  0.2R  -  26.12  ft 
Pitch  angle  at  0.2R  -  64.31  deg 
Density/2  =  p/2  ■  1  slug/ft3 


r/R 

dC  tf 

_ 1L 

dx 

CL 

tan  Si 

0.2 

0.122 

0.222 

1.982 

0.3 

0.310 

0.266 

1.305 

0.4 

0.533 

0.242 

0.963 

0.5 

0.732 

0.194 

0.754 

0.6 

0.850 

0.154 

0.612 

0.7 

0.876 

0.123 

0.519 

0.8 

0.782 

0.099 

0.432 

0.9 

0.559 

0.077 

0.372 

1.0 

0 

0 

0.324 

The  moments  of  thrust  and  torque  will  be  obtained  from  the  following  derivations  and 
by  the  use  of  Simpson's  Rule: 

For  the  moment  of  thrust  MT 


dT^-tfvVZdc 


Ti 


dT  dT  €*M08 

°  4  *  /». 


dTb  -^H2nVa2  Ul-e  tan  pjdC^ 
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Introducing  the  nondimensional  radius  x  -  r/R,  one  obtains: 

P 


dM 


Tb  -2R3nV°2l(1-  €  160  0<)(*“*o)dCrj 


MTh  -J-S3irV«:2  I  }  (x  -  *0)(1  -  «  tan  0,)  -Jj  rf» 


lb  2 


For  the  moment  of  torque  Mq 

dQ(-r  tan  (3 »  o?  2^ 

Introducing  the  nondimensional  radius  one  obtains; 


“ 2ff3"F°2  “  *o)  (tan  .  0,  +  e  )  dC^ 

(®-*0)(tan  f}i  +  ')<0!lidx 

b  2  *0 
The  constant  for-Af  and  A/p^  is 

/?3irV0>2  I  -  1,309,000  ft-Ib 


Further: 


f/ff 

c 

«  tan  0,- 

1  -  «  tan  0  { 

(tan  Pi+e) 

0.2 

0.036 

0.0714 

0.929 

2.018 

0.3 

0.030 

0.0392 

0.961 

1.335 

0.4 

0.033 

0.0318 

0.968 

0.996 

0.5 

0.041 

0.0309 

0.969 

0.795 

0.6 

0.052 

0.0318 

0.968 

0.664 

0.7 

0.065 

0.0337 

0.966 

0.584 

0.8 

0.081 

0.0350 

0.965 

0.513 

0.9 

0.104 

0.0387 

0.961 

0.476 

1.0 

bo 

- 

- 

- 
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The  following  calculations  are  restricted  to  the  hub  radius  *0 
radius,  follows  from  Simpson's  Rule: 


0.2.  Mr  ,  for  this 

V 


1 

2 

3 

4 

5 

r/R 

dCj 

(*-*o> 

Col. 

1x2 

Simpson's 

Multipliers 

f(MT^ 

0.2 

0.113 

0 

0 

1 

0 

0.3 

0.298 

0.1 

0.0298 

4 

0.1192 

0.4 

0.516 

0.2 

0.1032 

2 

0.2064 

0.5 

0.709 

0.3 

0.2127 

4 

0.8508 

0.6 

0.823 

0.4 

0.3292 

2 

0.6584 

0.7 

0.846 

0.5 

0.4230 

4 

1.6920 

0.8 

0.759 

0.6 

0.4554 

2 

0.9108 

0.9 

0.537 

0.7 

0.3759 

4 

1.5036 

1.0 

0 

0.8 

0 

1 

0 

2- 

=  5.9412 

M 


^09, .000  ^ 

3 


1,309,000  •  0.1  •  5.9412 
3 


259,230  ft-lb 


M0  is  obtained  by  the  same  method  as  applied  to  Mr 
V  5  b 


1,309,000  •  A*  •  2 
~~~  - 


164,950  ft-lb 


The  bending  moments  about  the  center  of  gravity  are  given  as  follows: 


-  Mr  cos  <f>  +  Mn ,  sia<f> 

*o  Tb 

sin  <f>  «  0.901 

cos  <f>  -  0,434 

Mx  -  259,230  •  0.434  +  164,950  •  0.901  -  261,130  ft-lb 
o 

-  Mr  sin  6  -  Mn  cos  <6 

y 0  Tb  *  Qb  * 

-  259,280  •  0.901  -  164,950  •  0.434  -  161,980  ft-lb 


The  stress  on  the  back  of  the  blade  at  point  of  maximum  thickness-is  given  by: 


*3 Myc 


y<> 


The  distances  from  the  center  of  gravity  to  the  point  of  maximum  thickness,  and  the 
moments  of  inertia  are  given  as  follows  (see  page  9) 


*3  « (o.47325  -  0.026 y  )l  -  0.43611  +  0.17  (y)  l 


0.184  ft 


y3  «  0.51  +  0.984/-  0.114  y+ 0.79  (-y)  l 

-  0.462  ft 

/  -  [b  ( 1/1 )2  +  0.04487  (t/l)3  ]14  b  from  Figure  5 

0 

-  [0.009(//1  )2  +  0.04487  (t/l  )31  lA 

-  0.1645  ft4 

I  -  (c  t/l)  lA  -  (0.04407  t/l) lA  c  from  Figure  6 
70 

-  4.958  ft4 

The  stress  on  the  back  at  the  point  of  maximum  thickness  is  thus: 

0.462  •  261,130  0.184  •  161,980 

0.1645  4.980 

-  -  739,400  psf 


or  5135  psi  compressive  stress. 

In  like  manner,  the  stresses  at  the  rose  and  tail  have  been  calculated  and  are  44.3  psi 
tension  and  1,179  psi  tension,  respectively. 
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For  complete  propeller  stress  analysis,  these  calculations  must  be  made  at  a  number 
of  radii. 
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